Structure of spanning trees on the two-dimensional Sierpinski gasket 
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Abstract 

Consider spanning trees on the two-dimensional Sierpinski gasket SG{n) where stage n is a non-negative 

Q_i, integer. For any given vertex x of SG{n), we derive rigorously the probability distribution of the degree 

' j G {1, 2, 3, 4} at the vertex and its value in the infinite n limit. Adding up such probabilities of all the 

vertices divided by the number of vertices, we obtain the average probability distribution of the degree j. 

The corresponding limiting distribution (pj gives the average probability that a vertex is connected by 1, 

2, 3 or 4 bond(s) among all the spanning tree configurations. They are rational numbers given as (pi = 

^ ■ 10957/40464, (/>2 = 6626035/13636368, (/>3 = 2943139/13636368, (/)4 = 124895/4545456. 
p 
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I. INTRODUCTION 



The enumeration of the number of spanning trees NsriG) on a graph G was first considered 
by Kirchhoff in the analysis of electric circuits ll^ . It is a problem of fundamental interest in 



28L 



3811 and physics [135 



The number of spanning trees corresponds to 



mathematics [|3L |4], 

a special g — limit of the partition function of the g-state Potts model in statistical mechanics 



a sp i 
111 



4111 . which in turn is related to the sandpile model 



131]. Just like other limits of the q- 



state Potts model, the spanning tree problem has been investigated intensely for decades, and 
has various applications in many areas. See, for example, ll39ll and references therein. It is also 
well known that there is a bijection between close-packed dimer coverings with spanning tree 
configurations on two related lattices ll36ll . Some studies on the enumeration of spanning trees 
and the calculation of their asymptotic growth constants on regular lattices were carried out in 



Refs. 115i|6, 



3711 . Once the total number of spanning trees and its asymptotic growth constant is 



obtained, the next step is to understand the geometric structure of spanning trees. One interesting 
question is the probability distribution of the degree of a certain vertex among all the spanning 
trees yj]. The geometric properties of spanning trees on Z*^ lattices, especially the square lattice, 
had been considered in 14, 3oll . 

Fractals are geometric structures of (generally non-integer) Hausdorff dimension realized by 



3. 



A 



repeated construction of an elementary shape on progressively smaller length scales yj 
well-known example of a fractal is the Sierpinski gasket that has been extensively studied in several 



contexts 



iwn example or a rractai is tne biei 

fl y y Q, y, y, Q, H 



23 



26 



27 



321]. Recently, the authors derived rigorously 



the number of spanning trees on the Sierpinski gasket and conjectured the result for arbitrary 



dimension Il7|]. It is of interest to consider geometric structure of spanning trees on self-similar 
fractal lattices which have scaling invariance rather than translational invariance. Different from 
the lattices that have translational invariance, e.g. the square lattice, the probability distribution of 
the degree on Sierpinski gasket depends on the vertex location. Thereby, it is natural to investigate 
the average of the probability distribution of the degree over all the vertices on SG{n) as n tends 
to infinite, and compare the values with the corresponding results on the infinite square lattice 
which is also 4-regular. In this paper, we shall present such probability distribution of the degree 
at any given vertex x on the two-dimensional Sierpinski gasket and the average, and the limiting 
distribution when the number of vertices goes to infinity. 
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II. PRELIMINARIES 



We first recall some relevant definitions for spanning trees and the Sierpinski gasket in this 
section. A connected grap h (without loops) G = {V, E) is defined by its vertex (site) and edge 



(bond) sets V and E |13i|22|]. Let v{G) = \V\he the number of vertices and e{G) = \E\ the number 
of edges in G. A spanning subgraph G' is a subgraph of G with the same vertex set V and an edge 
set E' C E. As a tree is a connected graph with no circuits, a spanning tree on G is a spanning 
subgraph of G that is a tree and hence e(G") = v{G) — 1. The degree or coordination number 
ki of a vertex Vi G V is the number of edges attached to it. A A; -regular graph is a graph with 
the property that each of its vertices has the same degree k. In general, one can associate an edge 
weight to each edge connecting adjacent vertices t>j and Vj (see, for example IsV]). For simplicity, 
all edge weights are set to one throughout this paper, so that the weight of each spanning tree is 
the same. 

The construction of the two-dimensional Sierpinski gasket SG{n) at stage n is shown in Fig. 
[B At stage n = 0, it is an equilateral triangle; while stage n + 1 is obtained by the juxtaposition of 
three n-stage structures. For the two-dimensional Sierpinski gasket SG{n), the numbers of edges 
and vertices are given by 

e{SG{n)) = 3"+^ , v{SG{n)) = ^(3" + 1) . 

Except the three outmost vertices which have degree two, all other vertices of SG{n) have degree 
four. In the large n limit, SG is 4-regular. 

Let us define the notation for the vertices of SG{n) to be used. An illustration for SG{A) is 
shown in Fig. |2l The denotation of the vertices is given progressively with increasing number of 
digits in the subscript as follows. First of all, fix o as the leftmost vertex. Consider the SG{m) with 
< m < n always has o as its leftmost vertex, and denote am and bm as its rightmost and topmost 
vertices, respectively. Cm is defined such that the vertices a„i, bm and demarcate the largest 
lacunary triangle of SG{m + 1). We then define the vertex in the middle of the line connecting 
am and am+i with m > 1 as am,i- Similarly, am,i and the associated bm,i and Cm,i demarcate 
a lacunary triangle with bm.i on the left and Cm.i on the right. Next for m > 2, we append the 
subscript m, 1, for the vertices of the largest lacunary inside the triangle with outmost vertices 
Om, ttmA, bm,i', the subscript m, 1, 1 for the vertices of the largest lacunary inside the triangle with 
outmost vertices a^,!, ctm+i, Cm,i', the subscript m, 1,2 for the vertices of the largest lacunary 
inside the triangle with outmost vertices bm,i, Cm,i, Cm, etc. In general for the vertices of SG{n), 



we use the notation where x = a,b,c and the subscript 7 = (71, 7s) has s components with 
l<s<r;,, l<7i<r2 and 7^ G {0,1,2} for k E {2,3,...,s}. For the vertices above the 
extended line connecting o and Cq, we will also use the notation x^ such that it is the reflection of 
the vertex with respect to this line. For examples, = ^21^ ^221 = ^212, C222 = £211, etc. The 
advantage of such vertex notation is that the quantities to be studied for the vertices ... with 
s > 2 components in the subscript can be expressed in terms of the quantities for the vertices with 
s — 1 components in the subscript as shown in Section IVl 



A ^ 

SG{0) SG{1) SG{2) 5G(3) 
FIG. 1: The first four stages n = 0, 1, 2, 3 of the two-dimensional Sierpinski gasket SG{n). 

Let us define the following quantities as in [9] . 

Definition II. 1 Consider the two-dimensional Sierpinski gasket SG{n) at stage n. (i) Define 
f{n) = NsTiSG{n)) as the number of spanning trees, (ii) Define g{n) as the number of spanning 
subgraphs with two trees such that the vertex hn belongs to one tree and the set of vertices {o, a„} 
belong to the other tree. (Hi) Define h{n) as the number of spanning subgraphs with three trees 
such that each of the outmost vertices 0, an and hn belongs to a different tree. 

Notice that for the spanning subgraph configurations counted by g{n), it is possible that the 
vertex 6„ is an isolated vertex with no bonds of trees connecting to it. A similar statement applies 
to the outmost verticex o, a„, 6„ for the spanning subgraph configurations counted by h{n). For a 
given vertex, we would like to investigate the number of bonds of spanning trees connecting to it 
among all the spanning tree configurations. We have the following definitions. 

Definition II.2 Consider the two-dimensional Sierpinski gasket SG{n) at stage n. For a certain 
vertex x G V{SG{n)), the number of bond(s) connects to it in a spanning tree configuration is 
denoted as j G {1, 2, 3, 4} or i G {0, 1, 2, 3, 4}. (i) Define fj{n, x) as the number of spanning 
trees such that there is (are) j bond(s) connecting the vertex x. Define the probability Fji^n, x) = 
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FIG. 2: The notation for the vertices of the Sierpinski gasket SG{A). The vertices inside the triangle (63, C3, 64) 
are reflection of the vertices inside the triangle (as, C3, 04) with respect to the line connecting o and C3, and are not 
shown. 

/,(n, x)/f{n). (ii) Define gi{n, x) as the number of spanning subgraphs with two trees such that 
the vertex bn belongs to one tree and the set of vertices {o, a„} belongs to the other tree, and there 
is (are) i bond(s) connecting the vertex x. Define the probability Gi{n, x) — gi{n, x)/ g{n). (Hi) 
Define hi {n, x) as the number of spanning subgraphs with three trees such that each of the outmost 
vertices o, an, bn belongs to a different tree, and there is (are) i bond(s) connecting the vertex x. 
Define the probability Hi{n, x) = hi{n, x) / h{n). 

For any vertex x of SG{n), the following relations for the probabilities should be satisfied, 

444 

Fj{n, x)^Y^ Gj{n, = ^ Hj{n, x) ^ 1 , 

j=l i=0 i=0 

which serves as a check for the results obtained. 

In this paper, we derive rigorously Fj{n, x) or fj{n, x) for an arbitrary vertex x e V{SG{ri)) 
with j — 1, 2, 3, 4. Such probability on translational invariance lattices in the infinite-vertex limit 
is independent of the vertex location. In contrast, as the Sierpinski gasket is a self-similar fractal 

5 



lattice which has scaling invariance rather than translational invariance, our results depend on the 
location of x. We shall consider the simplest vertex a; = o to obtain Fj(n, o) as Theorem IIII. 1 1 and 
its infinite n limit as Corollary IIII. II in Section [nil then move on to the vertices x E {am, bm,Cm} 
with < m < n to have Theorem IIV. 1 1 and Corollary IIV. 1 1 in Section HVl Fj{n,x) for the rest 
vertices will be treated in Section |V] as Propositions IV. II and IV.2I The summation and average of 
all Fj{n, x) for a given stage n will be studied in Section |Vll and such average in the infinite n 
limit will be obtained as Theorem IVI. II 



III. Fj(n,o) WITH J G {1,2} 



Consider the Sierpinski gasket SG{n) at stage n. We will derive Fj(n, x) for the vertex x = o 
in this section. Since the leftmost vertex o has degree two, fj{n,o) = for j = 3,4 and any 
n > 0. Similarly, we only need i E {0, 1, 2} for gi{n, x) and hi{n, x) with x E {o, a„, Due 
to the symmetry of SG{n), we have fj{n, o) = fj{n, an) = fjin, bn), gj{n, o) = gj{n, an) with 
j = 1,2, and hi(n,o) = hi{n,an) = hi{n,bn) with i = 0,1,2. According to the definition, 
go{n, o) = gain, an) = 0, but gain, bn) ^ for any n > 0. In fact, go{n, bn) is the only go{n, x) 
with non-zero value, and ho{n, x) is non-zero only when x E {o, an, bn}. The initial values for 
X = o at stage n = are /(O) = 3 with decompositions /i(0, o) = 2 and /2(0, o) = 1, g{0) = 1 
with decompositions gi(0,o) = 1 and (72(0,0) = 0, h{0) = 1 with decompositions /io(0,o) = 1 
and hj(0, o) = for j = 1, 2. We also have go{0, bo) = 1 and gj{0, bo) = for j = 1, 2. 

The following recursion relations was derived in [|7|] for n > 0, 



f{n + 1) = 6fin)^g{n) , 

g{n + l) = f{n)'h{n) + 7f{n)g{n)\ 

h{n + l) = 12 f (n) g (n) h{n) + Ug{n)^ 



(3.1) 



as illustrated in Figs. IMSI f{n), g{n), h{n) were solved exactly in [|7|] such that they satisfy the 
relation Sg^n)"^ = f{n)h{n). It follows that the second and third lines of (|3.1I) can be simplified as 



g{n + 1) = lQf{n)g{n)^ = ff{n)'h{n) , 
h{n + l) = b{)g{n)'^ = ^f{n)g{n)h{n) . 

Using Figs. [3l|5]for vertex o, we obtain the following recursion relations for j = 1,2: 



(3.2) 
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FIG. 3: Illustration for the expression of /(n + 1). 




FIG. 4: Illustration for the expression of g(n + 1). 



l,o,^4/,(„, ^^^^ 

^ 9j{n + 1,0) = fj{n, o)f{n)h{n) + 3/j(n, o)g{nY + Agj{n, o)f{n)g{n) , 

and 

gj{n + 1, bn+i) = hj{n, o)f{nf + fj{n, o)g{nf + 2[gj{n, o) + 2gj{n, hn)]f{n)g{n) , 
hj{n + l,o) = 4:hj{n,o)f{n)g{n) + 6gj{n,bn)g{nY + 4:fj{n,o)g{n)h{n) (3.4) 
+8gj{n, o)g{n)'^ + 2[gj{n, o) + gj{n, bn)]fin)h{n) . 

Setting /o(n, o) = and go{n, o) = 0, (|3.4I) reduces to 

5(o(n + 1, = ho{n, o)f{nf + Again, hn)f{n)g{n) , 
/io(n + 1, o) = 4/io(n, o)f{n)g{n) + 65fo(ra, hn)g{nf + 25(o(n, hn)f{n)h{n) . 



x3+ A---7.X3+ x3 + 



+ .<••■■>. x3+ .<•■■■>. x3+ .f--\x3+ /:•••■>. x3 



+ .<----A + 



FIG. 5: Illustration for the expression of h{n + 1). The multiplication for the eight configurations on the right-hand- 
side corresponds to three possible orientations. 
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(3.6) 



The initial values for the probabilities are 

Fi(0, o) = 2/3 , F2(0, o) = 1/3 , Gi(0, o) = 1 , 6*2(0, o) = 0, 
i7o(0, o) = 1 , Go(0, 60) = 1 , ^,(0, 60) = H,{{], o) = with j = 1, 2 . 

Divide the quantities in (|3.3I) - (|3.5I) by f{n + 1), g{n + 1) or h{n + 1) given in (|3.1I) or (|3.2I) . we 
get 



Fj(n + 1, o) = |F,(n, o) + \Gj{n, o) 
Gj{n + l,o) = |Fj(n,o) + |Gj(n,o) 



(3.7) 



Gj{n + 1, bn+i) = lGj{n, hn) + jQHj{n, o) + o) + \Gj{n, o) , 

+ 1,0) = 6„) + o) + |F,(n, o) + l^G,{n, o) 

for j = 1, 2, and 

G'o(^^ + l,&n+i) = |G'o(n,6„) + ^i7o(n,o) , 



(3.8) 



(3.9) 



i/o(n + l,o) = §Go(n,6„) + fi/o(n,o) . 

The probabilities Fj{n, o) and Gj{n, o) with j = 1,2 can be solved exactly by linear algebra as 
follows. 

Theorem III.l For the Sierpinski gasket SG{n) with non-negative integer n, 
H ) ; ^4 42^5^ ' ^ ' ^ 14 14^5^ ' 

2V ' ^ 14 42^5^ ' ' ^ 14 1445^ 

Proof Denote the vector Vj{n) = {Fj{n, o),Gj{n, o))^ with j = 1, 2 and 7i > 0. By (Il6l) 
and (|3.7I) . we have 

Vjin) = AW.iO) , 
where the matrix A and the initial Vj (0) are 





2 1 




2 




1 


A = 


3 3 


and V1(0) = 


3 




3 




3 2 
5 5 


1 






The matrix A can be diagonalized such that 

V,{n)=QADlQ~/V,{0) 
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where 



1 



and Qa 



1 5 
1 -9 



Therefore, 



V,{n) = QaDIQ-/V,{0) 



14 14^15' 14 14^15'' 
14 14^15'' 14 ~'~ 14^15'' 



for j = 1,2, and the proof is completed. □ 

From Theorem llll. 1 l and the exact expressions of f{n), g{n), h{n) in we have the following 
corollary. 

Corollary III.l For the Sierpinski gasket SG{n) with non-negative integer n, 

■11 5,1, 



fi{n,o) 



Ll4 42'15' 



f2{n,o 

probabilities for the vertex a are 



L14 ^ 42^15^ 



2o(n)2/3(n)g7(n) 



2«{n)2/3(n)g7(n) 



where a{n) = i(3" - 1), P{n) = 1(3"+^ + 2n + 1) and 7(n) = i(3" - 2n - 1). T/je //m/ring 



lim Fi(n,o) = ^ 



lim F2 (n, o) 



14 



In order to derive the probability Fj(n, x) for arbitrary vertex a; 7^ o, we need the following 
lemma: 



Lemma III.l For the Sierpinski gasket SG{n) with non-negative integer n, 



Go{n,bn) 


33 
~ 28 




^ 1 
^ ~ 28^25^ ' 






Ho{n,o) 


11 
14 




3 1 
' ^ 14^5^ ' 






Gi{n,bn) 


11 
~ 14 






14^ 


1) 

25^ 


G2in,bn) 


3 

~ 14 


+ 


9 1 Q 3 
7^5^ 28^^ 


5 

28 


^25 


Hi{n,o) 


11 
14 


+ 


3 1 4 3 
14^15^ ~7^5^ 


--( 


1) 

25^ 
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^ ' ^ 14 14^5^ 14^5^ 14^25' 



Proof Denote the vector Zo{n) = {Go{n,b„), Ho{n, o))"^ . By (|3.6I) and (13.91) . we have 
Zo(0) = (l,lf and 

Zo(n) = 5"Zo(n) = QbD],Q-'Zo{0) , 

where 

r „ „ 1 r „ 1 ^ 



Then 



B 



2 _3_ 

5 10 

_6_ _6_ 

25 25 



1 

IE 



Q 



B 



2 -6 



(3.10) 



Zn(n) 



14V5>' ~'~ 14^25^ 28 '^5'' 28 25 
7^5'' 7^25'' 14 V 5'' "'"14^25'' 



^o(O) 



28^5'' 28^25'' 
11 / 3\n I 3 / 1 \" 



14 V 5' 



14 V 25 ' 



gives Go{n, hn) and i/o(^, o). 

For the other probabilities, denote the vector Wj{n) = {Gj{n, bn), Hjin, o))^ with j = 1, 2, 



and split it into two parts Wj (n) = Wf^ (n) + Wf' (n) MM- By ^M, we have 



Wl^\n + 1) = 5iyf ^(n) and ivf ^(n + 1) = 5iyf V) + Rji^) , 



(3.11) 



where 



From Theorem lIII.il we know 



|F,(n,o) + ^G,(n,o) 



33 _|_ 13 / 1 

140 ~'~ 420 V 15'' 

143 I 11 / 1 

350 350 V 15'' 



and i?2 (?^) 



9 
140 

39 
350 



13 / 1 \n 
420 V 15'' 

11 / 1 
350 V 15'' 



To solve denote 



ai 



J_\n 
15) 



and W. 



(2), 



&i+&2(3^r 



15'' 



63 + &4( 



where the unknown Oj and 6j with i E {1, 2, 3, 4} satisfy the following relations by (13.111) : 



02 

15 

03 

a4_ 
15 



140 + s'^'l + 10^3 
420 + 5^2 + 10^4 

ii + + «3) 

+ ^(02 + 04) 



&1 

62 

15 



64 
15 



ifo + l&l + T7T&3 



10' 



13 I 2t I 3 I, 
~420 + 5''2 + 10"^ 



(3.12) 



11 _|_ _6_ 
350 ~'~ 25 



(62 + &4 
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It is straightforward to solve (13.121) and obtain 



14 7^15'' 
14 ' 14^15^ 



and W^in) 



_3_ _|_ 2(J_\n 

14 14 V 15 



(3.13) 



such that Wi'^\0) = (1/2, 1)^ and W^^\o) = (1/2, 0)^. 
As Wj{0) = {0,0f for j = 1,2, we have ^^^^^^(O) = 



-1/2, 



-1)^, <^(0) 



-1/2,0) 



T 



and W';^\n) can be solved by (ITTOl) . (I3ll1) as 



and 



(1), 



14 V 5'' "'"14^25'' 28^5'' 28^25'' 



I4V5'' 

3(3y 

7V5^ 



7V25'' 14V5>' 14V25^ 



'7^5'' 14V25>' 



4|'3^ 

'7V5> 



7V25' 



(3.14) 



14V25'' 28^5'' 28^25'' 
7^25) 14V5>' 14V25^ 



28 V 5 28 V 25 ^ 

14 V 5 7 ^ 14 V 25'' 



(3.15) 



Combining (|3.13I) - (I3.15I) . Gj{n, 6„) and Hj{n^ o) with j = 1, 2 are solved. □ 



IV. Fj{n + m + 1, Xn) WITH x G {a, b, c} AND n > 0, m > 

Consider the Sierpinski gasket SG{n + m + l) withn > 0,m > 0. We will derive Fj(ri + l,x„) 
with j E {1, 2, 3, 4} for the vertex s„ G {a„, 6„, c^} first, then Fj{n + m + 1, Xn) with arbitrary 
m > in this section. The corresponding ^-, (^ + 1, x„) and Hj{n + 1, x„) with j G {1, 2, 3, 4} will 
be used in the next section. Notice that Go{n + 1, Xn) = HQ{n + 1, x„) = for x„ G {a„, 6„, c„} 
as these vertices are not outmost vertices of SG{n + 1). 

For the Sierpinski gasket SG{n+l),we know fj{n + l,an) = /j (n + 1 , 6„ ) = /j (n + 1 , c„ ) and 
hj{n + 1, a„) = hj(n + 1, 6„) = hj(n + 1, c^) with j G {1, 2, 3, 4} because of rotation symmetry. 
From the definition of gj{n, x), we have gj{n + 1, 6„) = gj(n + 1, c„) but they are distinct from 
gj{n + 1, a„). Using Figs. [3ll5]for the vertex a„ or bn, we obtain the following recursion relations 

/i(n + l,a„) = 2fi{n,o)goin,bn)f{n), 

/2(n+ l,a„) = 2f2{n,o)goin,bn)fin) + 2fi{n, o)[gi{n, o) + giin,bn)]fin) 
+2fi{n,ofg{n) , 

< 

/3(n + l,a„) = 2/2(72, o)[5(i(n,o) +5(1 (ra,6„)]/(ra) +4/i(n,o)/2(n,o)5((n) 
+2/1 (n,o) [5(2 (n,o) + g2{n,bn)]fin) , 
^/4(n + l,a„) = 2/2(n,o)[5(2(n,o) +5(2(n,6„)]/(n) + 2/2(n,o)25((n) , 
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gi{n + l,an) = 2fi{n,o)go{n,bn)g{n) , 

g2{n + l,an) = fi{n,ofh{n) + gi{n,of f\n) + Afi{n,o)gi{n,o)g{ji) 

+'2[fi{n)gi{n,bn) + f2{n)go{n,hn)]g{n) , 
^3(n + l,a„) 2/2 (n, o) /i (n, o) /i(n) + 25(1 (n, 0)5(2 (n,o)/(n) 

+4[/2(n, o)c/i(n, o) + 0)0/2^ o)]«/(n) 

+2[/2(n, o)gi{n, 6„) + o)g2{n, hn)\g{n) , 
5(4(71 + 1, an) = /2(n, o)2/i(n) + /2(n, o) [25(2(71, 6„) + 45^2(71, o)\g{n) + 5(2(71, o)V(?^) , 



5i(7z+l,6„) = fi{n,o)ho{n,o)f{n) + 2fi{n,o)go{n,bn)g{n) 
+2gi{n,o)go{n,bn)f{n) , 

5(2(71 +1,6„) = [/2(7l,o)/lo("',o) +/l("',0)^l("',0)]/(?T') 

+2/1 (rz, 0)51 (rz, o)5(rz) + 51 (rz, o)V(^) 
+2[/2(rz, o)go{n, hn) + o)gi{n, hn)]g{n) 
+2 [5(2(77., 0)50 (t?^, &n) + gi{,n, 0)51(77., 6„)]/(7i) , 

< 53(71 +1,6„) = [/2(7Z,0)/7i(7Z,0) + /i(71,0)/72(7Z,0)]/(7z) 

+2[/2(7i, 0)^1(71, o) + /i(7i, 0)5/2(71, o)\g{n) 

+2[/2(7l, 0)5(1(71, 6„) + /i(7l, 0)5(2(71, 6n)]5'("') 
+2[5'2(7l, 0)5(1(71, hn) + 5'l("', 0)5(2(71, 6n)]/("') 

+25l(7^,o)52(7^,o)/(7^) , 
54(71 + l,6n) = f2{n,o)h2in,o)f{n) + 2g2{n,o)g2in,bn)fin) 

+2/2(71, o) [52(71, o) + 52(71, bn)]g{n) + 52(71, o)V(?^) , 



12 



and 



/ii(n + l,a„) = 4:fi{n,o)ho{n,o)g{n) + 2fi{n,o)go{n,bn)h{n) 

+Agi{n, o)ho{n, o)f{n) + 85-1 (n, o)gQ{n, bn)g{n) , 

h2{n + l,an) = 4:[f2{n,o)hQ{n,o) + fi{n,o)hi{n,o)]g{n) 

+2[/i(n, o)gi{n, 6„) + /2(n, o)go{n, bn)]h{n) 
+^[92{n, o)ho{n, o) + gi{n, o)hi{n, o)]f{n) 
+2/1 (n, 0)5-1 (n, o)h{n) + 6gi{n, ofg{n) 
+8[fi'2(n, o)goin, bn) + fi'i(n, o)gi{n, 6„)]fi'(n) , 

/i3(n + l,a„) = 4[f2{n,o)hi{n,o) + fi{n,o)h2{n,o)]g{n) 

+2[/2(n, 0)^1 (n, 6„) + o)5f2(n, &n)]/iN 
+4[gr2(n, o)/ii(n, o) + o)/i2(n, o)]/(n) 
+2[/2(n, o)gi{n, o) + /i(n, o)g2{n, o)]h{n) 
+8[5'2(n, 0)5(1 (n, 6„) + gi{n, 0)5(2(71, 6n)]fi'H 
+1251 (n, 0)52(71, o)5((n) , 

h4{n + l,an) = 4:f2{n,o)h2{n,o)g{n) + 2/2(71, o) [52(77, &n) +52(77,0)1/7(77) 

+452(77, 0)/72(77, 0)/(77) + 852(77, 0)52(77, 6^)5(77) 

+65(2(77,0)^(77) . 



Using the identity 85(77)^ = f{n)h{n), it follows that 



Fi(77 + l,a„) 

^2(77+ l,a„) 
^3(77+ l,a„) 
^ ^4(77 + 1, a„) 



■Fi(n,o)Qo(n,bn) 
3 

F2(n,o)Go(n,fc n j 

3 

F2(:n,u)[(;i(ii,u)-(;i(iijjn)] 
3 

F2 (n,o) [G2 (n,o)+G2 (n,b„)] 
3 



Fi(n,o)^ I Fi (n,o)[Gi (n,o)+Gi (n,6„)] 
3 ~'~ 3 



2Fi(ii,u)F2(n,o) I Fi(n,o)i(;2(»,o)-G'2(it,bn)] 



(4.1) 



+ 



3 

F2(n,o)^ 



+ 



Gi(77+ l,a„) 
G2{n + l,a„) 
^3(77 + I, an) 



Fi{n,o)Go{n,bn) 
5 ' 

3Fi(n,o)2 Fi(n,o)[2Gi(n,o)+Gi(n,6n)] , F2(n,o)Go(n,6n) , Gi(n,o)^ 

"T r; ^ 10 ' 



10 

3F2(n,o)Fi(n,o) , F2 (n,o) [2Gi (n,o)+Gi (n,6„ )] , G2(n,o)Gi(n,o) 
5 5 5 



(4.2) 



Fi(n,o)[2G2(«,o)+G2(n,bn)] 



r-n -L 1 n \ - 3-^2 (n,o)^ , F2(n,o)[2G2(n,o)+G2(n,&n)] , G2(n,o)2 
^ U-4(^77 -h 1, a„J ~ 10 ."^ in ' 



10 
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Gi{n + l,bn) 

G2{n + l,bn) 

G3{n + l,bn) 
GAn + l,bn) 



3Fi{n,o)Ho{n,o) , [Fi{n,o)+Gi{n,o)]Go{n,bn) 

10 5 ' 

3[F2 {n,o)Ho{n,o)+Fi (n,o) Jf i (n,o)] , Gi{n,of , Fi{n,o)[Gi(n,o)+Gi(n,b„)] 

10 10 5 

I F2(n,o)Go(n,bn) . Go(n,bri)G2in,o)+Gi{n,o)Gi{n,b„) 

3[F2 {n,o)Hi (ra,o)+Fi {n,o)H2 (n,o)] , [Gi (n,&n)+Gi (ra,o)]G2 (n,o)+Gi (n,o)G2 (n.fen) 

10 5 
I Fr(n,o) [G3-r(n,o)+G3-r{n,b„)] 

3F2{n,o)H2{n,o) ^ F2 (n,o) [G2 (n,o)+G2 (n,fe„)]+G2 (n,o)G2 (n,bn) _^ G2(n,o)^ 



(4.3) 



10 



10 



and 



Hi{n + I, an) 
H2{n + l,a„) 



3Fi{n,o)[2Ho{n,o)-\-Go{n,b„)] , 6Gi{n,o)Ho{n,o) , 4Gi(n,o)Go(n,b„) 

~'~ 25 ~'~ 25 ' 



25 



25 



6[F2(n,o)go(n,o)+Fi(n,o)gi(n,o)] 3Fi(ra,o)[Gi(ra,o)+Gi(n,fcn)] 

25 ~'~ 25 

. 3F2(n,o)Go(n,b„) 4[G2(n,o)Go(n,b„)+Gi(n,o)Gi(n,b„)] 

25 ~'~ 25 

I 3Gi(re,o)^ I 6[-go(n,o)G2(Ti,o)+gi(n,o)Gi(n,o)] 

25 ' 



25 



6[F2 (ra,o)Jj-i (ra,o)+Fi (n,o)ij"2 (n,o)] , 3F2 (ra,o) [Gi (ra,o)+Gi (ra,6n)] 
25 ~'~ 25 

3Fi(n,o)[G2(n,o)+G2(n,b„)] 6[G2(n,o)Hi(n,o)+Gi(n,o)H2(n,o)] 
25 ~'~ 25 

I [4Gi(n,b„)+6Gi(n,o)]G2(n,o) 4Gi(n,o)G2(n,fen) 

25 ' 



25 



Hi{n + I, an) 



6F2{n,o)H2{n,o) . 3F2 (n,o) [G2 («,o)+G2 {n,b„)] . 6G2 (n,o)H2 (n,o) 
~r oc; ~r 



25 ' 25 

4G2(n,o)G2(n,fen) , 3G2(rt,o)^ 



25 



+ 



(4.4) 



25 ' 25 ■ 

Next consider the Sierpinski gasket SG{n + m + 1) with n > and m > 0. The left-hand- sides 
of Figs. [3ll5]now represent SG{n + m + 1) with positive integer m, such that x„ G {a„, 6„, c„} 
locates within the lower- left triangle representing SG{n+m) in the right-hand-sides of the figures. 
As the vertices are denoted such that x^^ and are reflection of each other with respect to the 
extended line connecting o and cq, we have a„ = bn, bn = an and c„ = c„. We obtain the 
following recursion relations for j E {1, 2, 3, 4}: 

4/j(n + m, Xn)f{n + m)g{n + m) 



F,-(n + m + 



6/(n + m)'^g{n + m) 
^ [gfj(n + m, Xn) + gjjn + m, x„)]/(n + m)^ 

6/(n + m)'^g{n + m) 
2Fj(n + m, x„) (^^(n + m, Xn) + Gj{n + m, Xn) 

3 6 



(4.5) 



and 



Gj(?7, + m + 1, x„) = — — — h ' 



(4.6) 



5 10 10 

By symmetry, we know Fj{n + m, Xn) = Fj{n + m, x^) with G {a^, 6n, Cn} and (^^(n + 
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m, Cn) = Gj{n + m, c„) for positive integer m. Let us define the 3x3 matrix 

Fj{n + m + 1, an) Gj{n + m + 1, a„) ^^(n + m + 1, bn] 
Bj{n + m + l,n) = ir (n + m + 1, 6„) ^^(n + m + 1, 6„) ^^(n + m + 1, a„) (4.7) 

Fj(n + m + l,c„) Gj(r;, + m + 1, c„) (^^(n + m + 1, 
for non-negative integer n, m and j E {1, 2, 3, 4}. For m = 0, Bj{n + 1, n) has been obtained 
with elements given in (|4TT1)-(l431). By (|431 ) and (|431) . we have 



Bj{n + ni + l,n) = Bj{n + m, n)L' , 



for any m > 1, where 



L' 



2 3 3 

3 5 5 

1 A J_ 

6 10 10 

1 J_ _3_ 

6 10 10 



(4.8) 



We arrive at 



Bj{n + m + l,n) = Bj{n + 1, n)U'^ for all m > 0, n > . 
Solving Sj(n + m + 1, n) as in the proof of Theorem lIII.il its first column gives 

Fj{n + m + 1, a„) 



(4.9) 



Fj{n + m + 1, 6„) 
Fj{n + m + 1, c„) 



Q ^ 1 
^14 14^5^ ^ 



--ii-(-) 



+ l,a„) 
F,(n + l,c„) 



+ 1, a„) + (^^(n + 1, 6„) 
^^(ra + 1, 6„) + ^^(n + 1, a„) 
2Gj(n + l,c„) 



and we have the following theorem using (|4.1|) - (|4.3I) . 

Theorem IV.l For the Sierpinski gasket SG{n + m + 1) with non-negative integer n and m, 

P r-n -U m -U 1 n \ — I^IS ( _ _99_ ( i 55 ( 1 \^ 

riyiL -r III -r 1, an) — 5433 y^) ^372 1,25/ 16464 \375) 



(4)' 



121 



1 \ 



^3(72 + m + 1, a„ 



55 

5488 V5> 



('3^™ _ _22_ ('J_^™ _|_ 185 ( 1 y 

5488 V5/ 1029 V 25 / ~'~ 49392 V 375 / 

F („ 4- m 4- 1 n \ = 1^ — -§25. ('3\'" , rn_ ('J_^™ _ _55_ ( J_\"^ _ 2^ ('J_'|'" _L J_ ( J_'\ 
r2{ii^ '/t -r -L,Un; igg 548g 1^5^ "T -,^372 {25/ 5488 V 375 / 1176 V 15 / ~'~ 294 ^225/ 

i\m I 38 (±,\"^ _ 185 ('J_A™ _ 143 I 11 (J_\"^\ 

5/ ~'~ 1029 V25/ 16464 V375/ 3528 V 15 / ~'~ 2646 V225/ J ' 

^ 33 _ .855_ (3\"^ _ ^5_ ('J_^™ 1 _55_ (' J_^™ 1 H MV™ ]_ (' J-^l™ 

98 5488 V5/ 1372 V 25 / ~'~ 5488 V375/ ~'~ 147 V 15 / 147 V225/ 

(3\"^ _ 10 (±,\"^ _|_ 185 (J_\"^ _|_ M\™ _ 11 (' J_'|™\ 
V5/ 1029 V25/ ~'~ 16464 V375/ ~'~ 441 V 15 / 1323 V225/ J ' 

F (ri 4- m 4- 1 n \ = — ( ^'\'^ — -2^ H-Y" — 55 (' J_'|™ 1 JJ_ (l.Y' -^ J- /' J-) 

X'4i,/t-r //(. -r -L,Un; xgg 5488 1372 1,25/ 16464 ^75/ ~'~ 392 Vl5/ 294 V225/ 



57 
5488 



+ (4)' 



9 

5488 



(!)' 



2 

343 



('J.'l"^ _ 185 ('J_A™ _|_ JA /'J-'l™ -I- JJ_ rJ-^"^! 
U5/ 49392 U75/ 1176 Vl5/ 2646 ^225/ J 
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W (n ^ m ^^ r \ — - 22l. _L _5_ ( 

J-ll^/i -I- -t- 1, — 2744 1^5^ 343 V 25.' """ 8232 1375^ 

-L (l-Y^ [ ('3^™ _ 121 { _|_ 20 / 1 

~'~ V15/ \ 1372 V5/ 4116 V25/ ~'~ 3087 V375y J ' 

F in ^ \ r^-i21_i95_ /SN™ , _38_ M 5_ / J_^'" _ ^ M ^ 

-r2l,/t -r //t -r -L, 2744 V 5/ "^ 343 ^25/ 2744 V375/ 588 U5/ 

_|_ (IX^ /_55_ ('3\™ I _209_ (IV^ _ ('J_^™ _ 2!L {l-X^ _|_ 

~'~ V15/ \l372 V5/ 4116 V25/ 1029 V 375/ 441 V 15/ ~'~ 

Pr-n-Urri-Ul r \ — — -513- /'3\™- 10_ MV*" j 5_ /'J_^™ 1 Ji). /'J_^ 

j-sl^/i ±, C„; — gg 2744 343 V25;' "^ 2744 ^375;' "^147Vl5j 



1323 V225/ 



~'~ V15/ 1 1372 U/ 4116 V25/ ~'~ 1029 V375/ ~'~ 441 U5 / 1323 ^225/ 



, m 



pr-n-UrTi-Ul r \ — — (lY^ 6_ M ^™ 5_ ( J_Y^ 1 _5_ M V 

-r4i^/t -r //t -r -L, CnJ 2744 Uy* 343 ^25/ 8232 U75 """196 115^ 

_|_ il-Y'' ('3'|™' — JJ_ ('J-'l™' — i-^W^ -I- ^ ('J-'l™ -I- / 1 V 

~'~ V15/ \ 1372 V5/ 1372 V25/ 3087 V375/ ~'~ 147 V15/ ~'~ 1323 V225/ 

Corollary IV.l For the Sierpinski gasket SG{n + m + 1) with non-negative integer n and m, the 
limiting probabilities are 

lim„__oo + m + 1, x„) = , lim^^oo F2{n + m + l,Xn) = ^ , 

lim^^oo F'iin + m + 1, x„) = § , lim^^oo Fi{n + m + 1, x„) = , 

where the vertex x„ can be either an, or Cn. 

It is intriguing to notice that in Theorem lIV.il Fj{n + m + 1, a„) are distinct from Fj (n + m + 
1, Cn) with j E {1,2,3, 4}, while they have the same value in the infinite m limit. 



V. Fj{n + m, x) FOR GENERAL G V{SG{n)) WITH n > 0, m > 

Consider the Sierpinski gasket SG{n + m) with n > 0, m > 0. We will derive in this section 
Fj{n + m,x) with j E {1,2,3,4} for the general vertex E V{SG{n)) that has not been 
considered in previous sections. For the vertices inside the triangle with outmost vertices ctn-i, cin 
and Cn^i, let us append subscripts in the notation such that 7„ = (71 = n — 1, 72, ■ ■ ■ , 7^) with 
1 < s < n and 'jk £ {0, 1, 2} for k E {2, 3, s}. The results obtained in section HVl correspond 
to the vertex with s = 1 and n > 1, and we will tackle the vertex with s > 1 here. Similar to the 
definition of the vertex x^„^^, let us define the vertex as the reflection of x^^^ with respect to 
the line connecting a„ and 6„_i. By definition, we have 

where x can be either a, b or c, and 

Fj {n, ) = Fj (n, ) = Fj {n, %„^, ) , Hj (n, x^„_, ) = Hj {n, ) = Hj {n, ) , 
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due to the symmetry of SG{n). For m > 0, define the 3x5 matrix 



Fj{n + m,b^^J Gj{n + m,b^^J Gj{n + m,b^^J Gj{n + m,b^^J Hj{n + m,b^„J 
Fj{n + m,c^^J Gj{n + m,c^^J Gj{n + m,c^^J Gj{n + m,c^^J Hj{n + m,c^^J 



This is a generalization of Bj{n + m, n) in (14.71) . which corresponds to the case with s = 1. By an 
argument similar to that of ( 14.91 ). we have 



B,j{n + m, -fn,s) = Bj{n, -in,s)L' 



(5.1) 



for m > 0, where the 5x5 matrix 



2 3 3 3 _6_ 

3 5 5 5 25 

1 A X ^ 

6 10 10 50 

1 _!_ _3_ J_ J7_ 

6 10 10 10 50 

A A 

I 



is the generalization of L' in (14.81 ). It follows that the determination of Bj{n, -^n^s) for s > 1 will 
be sufficient. 

Let us first consider the vertices with s = 2 and 72 = 1, namely, 7„ 2 = (^—1,1) with 
n = 2, 3, .... We obtain the following equations (cf. Figs. | 



Fj(n,X(„_i,i)) 

Gj{n,X(^ri~i,i)) 
Gj{n,X(n~i^i)) 

Gj{n,xi^n~i,i)) 



2Fj{n-l,Xn-2) _|_ Gj{n-l,x„-2) , Gj{n—l,Xn-2) 
3 ~'~ 6 ~'~ 6 ' 



3Fj{n-l,x„-2) I 3Gj(n-l,Xn_2) , Gj(n-l,x„_2) 
5 10 ~'~ 10 ' 



3Hj{n-l,x„-2) _|_ Fj{n-l,x„-2) , 2G-, ( n-l,x„-2) 
10 ~'~ 10 ~'~ 



10 



Gj{n-l,x„-2) _|_ Gj(n-l,£n-2) 
10 

-l,Xn-2) I 3Gj(n-: 
5 10 '10 



3Fj{n-l,Xn-2) I 3Gj(n-l,£n-2) , Gj (n-l,a:n-2) 

5 10 
6Fj{n~l,x„-2) I 7Gj(n— l,a:n-2) , 6Gj(n—l,x„-2) 



25 



50 



+ 



25 



, 7Gj{n—l,x„-2) _|_ &Hj{n-l,x„-2) 



50 



25 



(5.2) 
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where x can be either a, b or c. Define the 5x5 matrix 



R 



2 3 

3 5 

1 3 



10 



25 



117 



10 10 10 50 







1 1 



2 
5 

1 3 



4 



6_ 

25 

7 



10 10 10 50 



0^0 



_6 
25 J 



(5.3) 



then (15.21) is equivalent to 



Bj{n, (n - 1, 1)) = Bj{n -l,n-2)R. 



(5.4) 



For general m > 0, we have the following formula combining (15.11) and (15.41) : 



Bj{n + m,{n- 1, 1)) = Bj (n, (n - 1, = Bj{n-l,n- 2)RL"' . (5.5) 



As Fj{n + m, i) = Fj{n + m, a;„_i,i) for x = a,b, c, the first column of the matrix in (15.5 
gives all Fj{n + m, j) in terms of the quantities for Xn-2- 

Proposition V.l For the Sierpinski gasket SG{n + m) with n > 2, m > 0, 



Fj{n + m,Cn-i,i) 



where 



Fj{n + m,dn-i,i) 
Fj{n + m, 6„-i,i) 
Fj{n + m,Cn-i,i) 

ei = (1,0,0,0,0)^ 



Bj{n-l,n-2)RL"'ei 



and] e {1,2,3,4}. 

Move on to the general vertex x^,^^^ ^ inside the triangle with outmost vertices a„, a„+i and 
c„ for the Sierpinski gasket SG{n + 1), where 7„+i,s = {n, 1,73, ...,7s) with 7^ e {0, 1,2}, 
k = 3,4, ..,s and 3 < s < n + 1. As 73 can take three possible values, let us discuss them 
separately. 

First consider the case with 73 = 1. The vertex ^ = a:(„ 1 1 -y^^.. .^^) is located inside the 
triangle with outmost vertices a„ 1, a„ + 1 and c„ 1. Associate with this , , . a vertex x^i 
where 7^ ^^-^ = (n — 1, 1, 74, 7^) has s — 1 component. That is, 7^ is obtained from ^n+i,s 
by taking out 73 = 1 and replacing 71 = n by n — 1. It can be seen that this vertex x^i _^ is located 
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inside the triangle with outmost vertices a„_i, a„ and c„_i. Moreover, x^i can be reached from 

in,s—l 

by a horizontal translation with the distance from a„ to o. Particularly, if s = 3 such that 



X. 



7n + l, 



7n+i,s = 1,1), then 7^ ^_i 



[n 



1, 1). By the method obtaining (15.41) , we have 



Bj (n + 1, 7„+i,,) = Bj (n, T^^^-i) R 



(5.6) 



if 73 = 1. 

Now consider the case with 73 = 2. The vertex x^^^^-^ ^ = X(^n,i,2,'Y4,...'ys) is located inside the 
triangle with outmost vertices 6„ 1, c„ 1 and c^. Associate with this x^.., , a vertex y^2 where 

^ ! IL-f-±.b in,s — l 

7n,s-i = ('^ ~ I5 1?74; ••■)7f)) has s — 1 componcnt. Here y = b when x = a and vice versa, 
and y = c when x = c. Namely, ?/ is related to x with three possibilities: (x, y) = (a, 6), (b, a) 
and (c, c). Similarly, 7^ is related to 7^ with three possibilities: (7^, 7^) = (1, 2), (2, 1) and (0, 0), 
where k = 4, s. Again, y;^2 ^ can be reached from x^^^^ ^ by a horizontal translation with the 
distance from a„ to o. We have 



BAn+ 1,7 



54n,7„^_i]i? 



(5.7) 



where 



Bj[n,^l,^^ 



Fj{n,b^2 ) Gj{n,b.2 ) Gj{n,b^2 ) Gj{n,b^2 ) Hj{n,b^2 ) 
F,(n, a^2 ) Gj(n,a^2 ) Gj(n,a^2 ) Gj{n,a^2 ) Hj{n,a^2 ) 

^J-(^'S^i.-l) GM^h.sj 



(5.8) 



By symmetry, the columns in (15.81) can be replaced as 



Fj{n,b^2 ] 
Fj(n, d^2 



Fj{n, c. 



in,s-l' 



Hj(n, a^2 

^ ' 'n,s-l ' 



Hj{n,b^2 ] 
Hj(n,a^2 

■J ^ ' 'n,s-l' 



and 



Gj{n,b^2 ) 
Gj{n,a^2 ] 
Gj{n,c^2 ] 

^ in.s — l'' 



Gj{n,b^2 ) 

^ in,s — l ' 

Gj(n,a^2 ) 

^ in,s — l ' 

Gj{n,c^2 ) 

^ in,s — l ' 



Gj{n,b^2 ) 
Gj{n, d^2 ] 
Gj{n,c^2 ] 

^ 'n.s — 1 ^ 



Gj{n,b.2 ) 

^ in,s — l ' 

Gj{n, d^2 ) 
Gj{n,c^2 ) 



so that 



^i(^>7n,s-l) = ^i(^>7n,s-l)^2 
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where 



1 

10 

10 

1 

1 



(15.71) can be rewritten as 



Bj{n + 1,7- 



n+l,s) 



if73 = 2. 

Finally consider the case with 73 = 0. The vertex x^^^-^ ^ 



^(n, 1,0, 74, .■■7s) 



(5.9) 



is located inside the 



triangle with outmost vertices a„, a„ 1 and 6„ 1. Associate with this x^^^-^ ^ a vertex ^^0 ^ where 
7° s-i — ('^ ~ I5 745 •••5 7s) has s — 1 component. Here z = c when x = b and vice versa, and 
z = a when x = a. Namely, z is related to x with three possibilities: (x, z) = (a, a), (6, c) and 
(c, 6). Similarly, 7° is related to 7^ with three possibilities: (7^,7°) = (1,0), (0, 1) and (2,2), 



4, s. We use the notation such that the vertex zm is the reflection of the vertex 

IJI.S—I 



where k 

^-y-^.s-i with respect to the line connecting a„_i and 6„. It can be seen that ^ can be reached 
7n+i s by a horizontal translation with the distance from a„ to o. We have 



from X 



Ej(n + l,7„+i,,) = 5j(ra,7°,_Ji? 



(5.10) 



where 



Fj(n,c^o 
Fj(n,6^o 



F,(n, a^o ) G',(n, a^o ) GAn^a^o ) Gjin^a^ 

) Gj(n,c^o ) Gj(n,c^o ) GAn.c^fi 

) Gj(n,b^o ) Gj(n,b;M ) Gj{n,b^) 

) Gj(n,a^o ) Gi(n,d;M ) Gj(n,a^o 

) Gj(n,c^o ) Gj(n,c^o ) Gjin.c^ 

) Gjin.b^o ) Gj(n,bM ) Gj(n,b^o 



Fj{n,b^ 



Hj (n, a^o 
Hj{n,c^o 

^ 'n.s — l 

Hj(n, 6^0 



Hj (n, a^o 

^ !n,s — l 

Hj(n, c^o 

^ in,s — l 

Hj(n, 6^0 

^ 'n,s—l 



by symmetry, so that 



^i(^,7n,s-l) = ^j(^,7n,s-l)^0 , 
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where 



En 



1 

10 

1 

10 

1 



(15.101) can be rewritten as 



Bj{n + 1, 7^+1,,) = Bj{n, lls-i)EoR 



(5.11) 



if73 = 0. 

Denote Ei = Jsxs as the identity matrix. (|5.6I) . (15.91) and (15.1 II) can be combined to give 

B,{n + l,^n+i,s)=B,{n,f,^^,_,)E,^R, (5.12) 

where 73 G {0, 1, 2}. As Fj{n + 1, x^^^^ J = Fj{n + 1, x^,^^-^ J for x = a,b, c, the first column 
of the matrix in (|5.12|) gives Fj{n + 1, a;^^^^ J for any vertex a;^„^^^ in terms of the quantities for 



X. 



7n,s-l • 



Proposition V.2 For the Sierpinski gasket SG{n + 1) with n > 2, consider the vertex a;^„+i 
where 7n+i,s = {n, 1, 73, 7^) with 3 < s < n + 1 and 7^ G {0, 1, 2} for k G {3, 4, .., s}. 



Bj{n,^l%_^)E^^Rei . 



Using Theorems IIII. 1 [ HvTI and Propositions IV. 1 1 |Y2] repeatedly. Fj{n + 1, a:^„+i,J for all the 
vertices of SG{n + 1) can be obtained. 



VI. SUMMATION AND AVERAGE OF Fj (n, x) OVER ALL THE VERTICES OF SG{n) 

It is worthwhile to derive the summation of Fj{n, x) over all the vertices x of SG{n), defined 



as 



x£V{SG(n)) 



21 



and the average of Fj (n, x) over all the vertices, defined as 



^'^""^ v{SG{n)) |(3" + 1)- 
It is clear that for any non-negative integer n, 

4 

For the vertex x^^ ^ with s = 1, i.e. am, bm and Cm, define their sum 

Xj{n, m) = Fj{n, am) + Fj{n, bm) + Fj{n, Cm) ■ 
Similarly for the vertex with s = 2, define 

Yj(n, m') = Fj{n, am',i) + Fj{n, bm',i) + Fj{n, Cm',i) , 
where m > 0, m' > 1 and n is larger than m and m'. By (15.11) . we have 

Xj{n,m) = {l,l,l)Bj{n,m)ei = (1, 1, l)5j(m + 1, m)L"-"'-^ei , 

y,(n, m') = (1, 1, l)Bj{n, (m', l))ei = (1, 1, l)5,(m' + 1, (m', 
= {l,l,l)Bj{m',m' -l)RU'~"''-^ei. 

The first few are 

$,(0) = 3F,(0,o), 

$,(1) = 3F,(l,o)+X,(l,0) = 3F,(l,o) + 2F,(l,ao) + F,(l,co), 

and 

$,(2) = 3F,(2,o)+X,(2,0)+X,(2,l) + 2yj(2,l) 

= 3Fj(2,o) + (1, 1, 1) |5, (1,0)L + 5,(2,1) + 25, (1, 0)i?}ei 

The corresponding values for j G {1,2,3,4} are 

$l(0) = ^, $2(0) = ^, $3(0) = $4(0) = , 

1 19 7 1 

$i(l) = -, $2(1) = ^ $(1) = ^ $4(1) = _ 

u ; 2 ' ^ 54 ' ^ ^ 54 ' ^ ^ 54 ' 
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^ ,o^ 163 ^ 5257 
*i(2) = 7^ , *2(2) 



450 ' ^ ' 12150 

For n > 3, we need the summation 



^ 2203 ^ 289 

$3(2) = , $4(2) 

^ 12150 ' ^ ^ 



12150 



n 2 2 



s=3 7s=0 73-1=0 73=0 



for the vertices with s > 3. By (|5.4I) and (|5.12l) . M, (3) = Xl73=o 73,3) is given by 

M,(3) = [b^{2,%)Eo + B^{2,%)E, + B^{2,%)E2}r 
= B,{2,{l,l))[Eo + Ei + E2]R 
= Bj{l, 0)R[Eo + Ei + E2]R = Bj{l, Q)RER , 



where 



E = Eq + El + E2 



3 

2 10 

1110 

12 

3 



The general expression for n > 3 is 

M,(n + 1) 

n+l 2 2 2 

s=3 7s=0 7s-i=0 73=0 

= [B,{n,t,2)Eo + B,{n,f^^,)Ei + i?,(n, 7^2)^2}^? 

n+l 1^ 2 2 2 



s=4 ^7s=0 7s-i=0 74 =0 

n 2 2 2 

[i?,(n,(n-l,l)) + X:{E E ••■Ei^.l 

s=3 7s=0 7s-i=0 73=0 



R 



ER 



\^Bj{n - l,n - 2)R + Mj{n)^ER 

Bjin - 1, n - 2)RER + Bj{n - 2, n - 3)R{ERf + Mj{n - l){ERf 



n-l 



(6.3) 



m=l 
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For example, 



$j(3) = 3F,(3,o) +X,(3,0) +Xj(3,l) +X,(3,2) 

+2[r,(3, 1) + rj(3, 2) + (1, 1, 1)5,(1, 0)RERei] 
= 3Fj{3,o) + (1,1,1)|5,(1,0)L2 + 5j(2,1)L + 5,(3,2) 

+2[Bj{l, 0)RL + Bj{2, 1)R + 5,(1, 0)RER]^ei 
= 3Fj(3,o) + {l,l,l)^^Bj{l,0)[L^ + 2RL + 2RER] + Bj{2,l)[L + 2R] + Bj{3,2)^ 

For general ri > 3, we have 



ei 



n-l 



n-1 



$,(n) = 3F,{n,o) + ^ Xj(n,m) + 2 ^ Yj{n,m) + 2(1, 1, 1)|^ Mj(m)L"-™}ei , (6.4) 

m=0 m=l m=3 

with Xj(n, m) and yj(n, m) given in (16.11 ) and (16. 2k respectively. From (16.31 ). the summation in 
the last term of (16.41 ) is 



X^M,(m)L- 



m=3 



n m— 2 
m=3 s=l 


[s,s-l)R{ER)"'- 


n— 1 m— 1 
m=2 s=l 


[s, s - l)R{ERr- 


n-2 

^5j(s,s- 


n-l 

-1) [r{er) 

m=s+l 


n-2 

^5^(s,s- 

s=l 


n—s—1 

- 1) j2 [Ri^^y 

m=l 



L 



n—l—m 



so that $j (n) can be calculated exactly for any positive integer n. 

Proposition VI. 1 For the Sierpinski gasket SG{n) with n > 3, the summation ofFj{n,x) over 
all the vertex is given by 

n n—1 



$j(n) = (1,1,1)<^ ^5j(m,m-l)i?"-'" + 2^5j(m,m-l)i?L 



n— m— 1 



m=l 



m=l 



n-2 



n—m—1 



m=l 



+2^5j(m,m-l) J2 iRiERYL''-^-"'-' [ei + 3Fj{n,o 



s=l 



(l,l,l)l^5,l 

^m=l 



m, m - l)i?"~" + 2Bj{n -l,n- 2)R 



n-2 



n—m—1 



+2^5j(m,m-l) J2 RiERyV-"^-^-' Ui + 3Fj(n,o). (6.5) 



m=l 



s=0 
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Let us consider limiting distribution 



lim (j)j{n) = 

n— >oo 



withj G {1,2,3,4}. It is easy to see that the term 3Fj(n, o) in (16.51) can be neglected in the infinite 
n limit for (pj, namely, 

lim ' ' = 



n^oo |(3" + 1) 

as the value Fj(n,o) is between and 1 for the four possible j. Similarly, the values of the 
quantities Fj{m, Xm-i), Gj{m, Xm-i), Hj{m, Xm~i) with x = a, b,c in the matrix Bj{m, m — 1) 
are between and 1, and all the eigenvalues of R given in (|5.3I) are positive and less than or equal 
to 1, such that 



lim 

n— >oo 



(1, 1, ^){EI=i B,im, m - 1)R--^ + 2B^{n - 1, n - 2)i?}ei 3^ + g 

< lim ^r— = 



K3" + 1) 



1(3" + 1) 



Therefore, only the double summation term in (16.51) gives non-zero contribution for (pj. Rewrite 
EE = Qi[Di + Di]Q^^ and L = Q2^2Q^\ where 



Qi 



159 -87 -3 -3 

38 14 1 2 1 

38 14 2 -4 

38 14 -1 2 1 

15 45 -3 5 



Q2 



18 -27 -2 
5 -1 98 -1 1 
5 -32 1 1 
1 -52 
13 



and 






















3 
















1 



















1 

































3 
10 



















2 

5 








, D,= 

























6 
25 



















3 

25 

































1 

5 



















0_ 

































1 

15 . 



then 



n—m—l 



n—m—1 



R{ERyL 



n—m—l—s 



R J2 [Qi{Dl + Dl)Qi' 

s=0 

R J2 \QiDIQi^]\Q2D 



n— m— 1— s^ — 1 
V2 



25 



n—m—1 



+R J2 [QiDlQi' 

s=0 

The elements of Bj(m, m — 1) have been solved in (I4.1I) - (I4.4I) . Define 

Zj{m) = (1, 1, l)Bj{m, m — 1) 



(6.6) 



(l)^3^ 
5^ 



25' 



+ A 



^375' 



for J = 1 



363 \ 1 \(2) 
196^0 + 


(1 


)™ + A^^^ 


He 


)" + A^ 


^375 


T + Af 


(1^ 


+ Af 


V225 


99 X 1 \(3)f 
98^0 + Ai [ 


1) 


" + Af ( 




™ + Af ( 


375'' 


™ + Af ( 


A) 


+ Af ( 




I 196^0 + Ai 


(1 


+ A?) 




+ Ai') 


V375 


r + 


(1^ 


r + A^') 


V225 



(6.7) 



where Aq = (1, 1, 1, 1, 1) and 

aS^ 



A 



(1 



,(2 



A 



(2 



(3: 



Ai' 



A 



(3: 



(3: 



A 



A? 



,(4 



605 121 121 121 


1089 




392' 56 ' 56 ' 56 


' 392 '* 




3125 -375 -375 


-375 


5 


392 ' 56 ' 56 


' 56 


'3 


2375 95 95 95 - 


-2109 




196 ' 28' 28' 28' 


196 




-1265 -187 -IS 


7 -18 


7 



),A(^) = ( 



-1375 -55 -55 -55 122L 



196 



(2) 



' 28 28 28 196 ' 
-275 -55 -55 -55 -495 



.(2) 



^ 396 ' 56 ' 56 ' 56 ' 392 ' 
-9375 1125 1125 1125 -1755, 



392 



196 ' 196 ' 196 ' 196 ' 196 



).Ar'=( 



56 ' 56 ' 56 
,10 -240 -240 -240 54 



392 



49' 49 ' 49 ' 49 '49^' 



-285 -57 -57 -57 -513 



),Ar = ( 



392 ' 56 ' 56 ' 56 ' 392 



-625 -25 -25 -25 555. 



-9375 -1125 -1125 -1125 1755 



392 ' 56 ' 56 ' 56 ' 392 " ' ^ 



196 ' 28 ' 28 ' 28 ' 196^ ' 
(3) _ 230 34 34 34 -162 

~ *^l9"'49'49'49' 49 ' 



500 480 480 480 -108, 



Al 



(4) 



49 ' 49 49 49 49 
-375 -15 -15 -15 333 (4) 

196 ' 'W' 'W' 196'* ' ^ 
345 51 51 51 -243 



-45 -9 -9 -9 



196 196 196 196 196 



A 



(4) 



392 ' 56 ' 56 ' 56 ' 392 ' ' 
-3125 375 375 375 -585. 

^ 392 ' ^' '5Q' 392 ' 
,250 -240 -240 -240 54, 



49 49 



49 



49 49' 



Substituting (1631) into (1631) . we get 

A "~2 

<Pj = Jim -(3" + l)-i 5^ Z,(m)i?{ J2 [QiDlQ^'][Q2Dr'^-'-V2' 

m=l s=0 

n—m—1 

+ J2 [QiDlQ^'][Q2Dr'^-'-V,']] 



n— m— 1 



(6.8) 



As the eigenvalues of Di and D2 are between and 1, the first term in (16.81) makes no contributions 
since 



n-2 



n—m—1 



m=l 

n-2 



s=0 

n—m—1 



ei 



m=l 



s=0 



ei < 3n . 
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Consider the second term in (16.81) . 

n— 2 n— m— 1 



m=l 
n-2 



s=0 



m— 1 



771=1 S=0 



Q2D\ 



Q2 



771=1 

n~2 



777=1 



where 



1 











and define 

D{n, m) 



n—m—l 
^ 3S^n-m-l- 

s=0 

'^s ( 6 ^n—TTi— 1— s 



71— m— 1 

E 

s=0 



25 [3" 










3^(ii 





1 Nn— »n— 1- 












3^ 









3S(i^77-r«-l- 



69 












isp^-^-Cj^)"-'"] 

44 
























5 [3 



14 











os( 3 Nn— 771— 1— ; 






10[3"-"'-(^)"-'"] 
27 



T-'^V + V2{n - m) 



27 



with 



V 



§ 

i 

I 

^0 

if 



m] 










_ ( 1 \n — m — 1 
IS-' 

44 










2 











14 











_ / _3_ \n — m—l 



Since the absolute values of all the eigenvalues of 1^2 — ^) are less than one, we have 



lim -(3" + i: 



n-2 



m=l 
-1, 



+Zj{m)RQiDiQ-'Q2V2in - m)Q 



1-1 



ei 



n-2 



lim -(3" + 1)-^ y 3"-'™Z,(m)^ei 



(6.9) 



m=l 



where R = RQiDiQ^^Q2'DQ2^ ■ Substituting the expression of Zj{m) from (16.71) into (16.91 ). 
carrying out the summation and taking the infinite n limit, we arrive at 



4 



(1) 



4 

363Ao 



2 I 3 
'W ~'~ 1124 
(2) ,(2l 



-Rei for j = 1 



\(2) -,{2) 

_|_ ±L l_ l_ ^^3 l_ ±i l_ ^5 

392 ~'~ 4 ~'~ 74 ~'~ 1124 44 ~'~ 674 

196 ~'~ 4 74 1124 44 674 



27Ao 



392 



+ 



A 



(4) 



+ 



Ai 



(4) 



74 



+ 



a: 



(4) 



1124 



+ 



A 



(4) 



44 



+ 



a: 



(4) 



674 



Rci for J = 2 
Rci for j = 3 , 
Rei for j = 4 . 



The matrix productions can be done to give the following theorem. 

Theorem VI.l Consider all the vertices of the Sierpinski gasket SG{n) in the infinite n limit. The 
average probabilities that a vertex is connected by 1, 2, 3 or 4 bond(s) among all the spanning 
tree configurations are 

10957 6626035 

= —— = 0.270783906682 ■ ■ ■ , 02 = = 0.485909077842 ■ • • , 

^ 40464 ' ^ 13636368 
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0.215830124267 ■ ■ ■ , 6^ = ——— = 0.0274768912073 ■ ■ 



13636368 ' ^ 4545456 

Corollary VI.l Consider all the vertices of the Sierpinski gasket SG{n) in the infinite n limit. 
Denote 9 as the average number of bonds connecting to a vertex among all the spanning tree 
configurations, then 

^ = 01 + 202 + 303 + 404 = 2 . 

We list the numerical values of 0j(n) with j G {1, 2, 3, 4} for < n < 5 and infinite n limit in 
Table m We find that 0i(n) decreases monotonically as n increases, while 02 (^), 03 ('^) and 04 (n) 
increase monotonically. The values for n = 5 are already very close to 0j in the infinite n limit 
with deviations about 1%. 

It is interesting to compare the Sierpinski gasket SG{n) in the infinite n limit with the infinite 
two-dimensional square lattice which is also a 4-regular lattice. For the square lattice, all the 
vertices are identical due to the translational invariant. The probabilities that a vertex is connected 
by 1, 2, 3 or 4 bond(s) among all the spanning tree configurations have been solved exactly in 
BOll that were denoted as fj with j e {1,2, 3, 4}. As shown in Table H /i, /s and /4 are slightly 
larger than 0i, 03, 04, respectively, while /2 is smaller than 02. Especially, the average number 
of bonds connecting to a vertex among all the spanning tree configurations on the square lattice 
fi + 2/2 + 3/3 + 4/4 is equal to two, which is exactly the same as 9 here for the Sierpinski gasket. 
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